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We apply time-dependent R-matrix theory to study inner-shell ionization of C atoms in ultra-short 
high-frequency light fields with a photon energy between 170 and 245 eV. At an intensity of 10^^ 
W/cm^, ionization is dominated by single-photon emission of a 21 electron, with two-photon emission 
of a Is electron accounting for about 2-3% of all emission processes, and two-photon emission of 
21 contributing about 0.5-1%. Three-photon emission of a Is electron is estimated to contribute 
about 0.01-0.03%. Around a photon energy of 225 eV, two-photon emission of a Is electron, leaving 
in either Is2s2p® or ls2p‘* is resonantly enhanced by intermediate ls2s^2p® states. The results 
demonstrate the capability of time-dependent R-matrix theory to describe inner-shell ionization 
processes including rearrangement of the outer electrons. 
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I. INTRODUCTION 

Over the last decade, great strides have been made in 
the development of free-electron lasers operating in the 
VUV - X-ray regime. Several free-electron lasers operat¬ 
ing in the XUV-X-ray regime have become available to 
the community in recenWears: for example, FLASH [l|, 
LCLS 0, and SACLA Q. These facilities have demon¬ 
strated their potential for opening new areas of atomic, 
molecular and optical physics, for example through the 
study of Auger resonances which cannot be excited by 
a single photon Q, multiphoton sequential ionization of 
Xe up to Xe^®"*" at a photon energy of 1.5 keV Q and 
multiphoton multiple ionization of N 2 

Photoionization in high-frequency laser fields tends to 
be dominated by the innermost electron that can be 
ejected. However, the outer electrons will also experi¬ 
ence the light field, and can therefore still absorb a pho¬ 
ton. Hence, a full description of the atomic or molecular 
response should consider all electrons that could possibly 
be affected by the laser field. In addition, outer electrons 
do not necessarily remain in their original orbital when 
an inner electron is removed from the system. The po¬ 
tential seen by the outer electrons may change suddenly, 
leading to shake-up excitation of the outer electrons. 

A full theoretical or computational study of the inter¬ 
action between high-frequency laser light and atoms will 
therefore require a method which can describe the simul¬ 
taneous response of many electrons to the laser field. Sev¬ 
eral such methods have been developed in recent years, 
such as the time-dependent configuration-interaction sin¬ 
gles (TDCIS) method, which has recently been applied 
to study above-threshold ionization for light elements in 
the hard X-ray regime Q, a Green’s function technique, 
algebraic diagrammatic construction (ADC), which has 
been applied to study fast dynamics in glycine using 
laser pulses at a photon energy of 275 eV and time- 
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dependent R-matrix theory, which has been used to study 
the competition between emission of a 2s and a 2p elec¬ 
tron in C in the UV regime 

In the present study, we continue our study of the re¬ 
sponse of C atoms to laser light by investigating the pho¬ 
ton energy range between 170 and 245 eV. The response 
of the carbon atom is of particular interest as it is the 
prime constituent of biological molecules. The removal 
of inner electrons from a carbon atom can provide new 
insight into molecular systems. For example, in [13, it 
was demonstrated that the removal of a Is electron from 
the carbon atom in methane could be exploited experi¬ 
mentally to extract information about the molecular ge¬ 
ometry. In 13, it was also proposed that dynamics in 
ionized glycine could be studied in a pump-probe scheme 
where the probe pulse excites a localised Is electron of 
C to the orbital in which a hole is created by the pump 
pulse. A photon energy of 275-280 eV was suggested for 
this purpose. Advances in laser technology have very 
recently been exploited to generate individual subfem¬ 
tosecond pulses in this photon energy range m- Thus it 
is of interest to investigate ultra-fast dynamics involving 
inner-shell electrons. 

As a first step towards the treatment of short laser 
pulses at a photon energy of 284 eV m , we compare in 
the present study emission of the inner Is electron with 
emission of the outer 2s or 2p electrons for the photon 
energy range between 170 and 245 eV. Numerical studies 
of ultra-fast dynamics at 284 eV involving the C atom 
will require great care with the pulse shape to ensure 
that inner-shell ionization processes are not dominated 
by direct single-photon emission arising from the outer 
edges of the pulse bandwidth. In the present photon en¬ 
ergy range, emission of a Is electron requires absorption 
of (at least) two photons, whereas the emission of a 2s or 
2p electron requires absorption of a single photon only. 
This comparison is similar to a previous comparison of 
two-photon emission of the Is electron versus one- and 
two-photon emission of the outer 2s electron using the 
R-matrix Flog net approach for Li“ [l3 or the Is2s 
state in He |13l |. 
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To study the response of the carbon atom, we use the 
recently developed time-dependent R-matrix approach 
RMT [3 t 13- If combines the capability of R-matrix 
theory to describe a wide range of processes in general 
atomic systems fT7l - [l^ with the computational capabil¬ 
ity of the HELIUM approach [ 2 ^. The combination of 
these two techniques has enabled the determination of 
time delays in Ne llol. and high-harmonic generation at 
mid-IR wavelengths [^. For these studies, relatively lit¬ 
tle atomic structure was taken into account. By adopt¬ 
ing an R-matrix with pseudo-states philosophy or 
an intermediate-energy R-matrix approach [23| , we have 
recently demonstrated that the RMT codes can also be 
used for the study of double photoionization processes 
To describe the double continuum accurately, 
extensive atomic structure needs to be taken into ac¬ 
count. The success of these latter studies suggests that 
the RMT approach is capable of treating atoms in strong 
fields with a detailed description of atomic structure. 

In the present study, we wish to explore the applica¬ 
tion of RMT theory to a case where electrons can be 
ejected from different shells with significantly different 
binding energy. This application poses new demands on 
the computational accuracy. The continuum needs to 
be accurate up to very high energies to describe rele¬ 
vant above-threshold ionization processes involving outer 
electrons. This has the computational consequence of in¬ 
creasing the round-off error in matching the wavefunc- 
tion at the boundary between the inner and outer R- 
matrix regions, where the description of the ejected elec¬ 
tron changes from basis-set techniques to finite-difference 
methods. Hence the application of the RMT approach to 
inner-shell processes presents new demands on the com¬ 
puter codes. 

In section m we give a short overview of the RMT 
approach. We also provide a brief description of the basis 
set used to describe the C atom, and the laser pulse. The 
results are presented in section uni with an emphasis on 
the competition between the multiphoton emission of an 
inner Is electron and single- and multiphoton emission of 
an outer 2£ electron. 


II. COMPUTATIONAL METHODS 

Time-dependent R-matrix theory is the extension of 
the standard R-matrix approach to the solution of 
the time-dependent Schrddinger equation [l3 - [l8l 
Although the initial applications of time-dependent R- 
matrix theory described electrons restricted to a finite 
region surrounding the nucleus [i3[il, subsequent im¬ 
plementations adopted the standard R-matrix concept 
of division of space into two distinct regions: an inner 
region in which all electrons are close to the nucleus, 
and an outer region in which one electron has moved 
well away from the others. In the original formulation 
of time-dependent R-matrix theory, an R-matrix propa¬ 
gation scheme was employed to propagate the wavefunc- 


tion [ 2 ^ . This approach relies on the solution of systems 
of equations throughout the calculation, and as a conse¬ 
quence calculation time increases rapidly with an increase 
in atomic structure. 

The most recent implementation of time-dependent R- 
matrix theory is R-matrix theory with time dependence 
(RMT) In this approach, the wavefunction in the 

inner region is described in terms of a standard R-matrix 
basis with a B-spline representation of the continuum or¬ 
bitals. The wavefunction in the outer region is described 
in terms of a direct product of a residual-ion state cou¬ 
pled with a finite-difference representation of the wave- 
function for the outer electron. Near the boundary be¬ 
tween inner and outer region, the wavefunction must be 
shared by the inner and outer regions. This is achieved 
through evaluation of the inner-region wavefunction on 
an outer-region grid extension into the inner region for 
use by the outer region. The outer-region wavefunction 
information needed by the inner region consists of spa¬ 
tial derivatives of the outer-region wavefunction at the 
inner-region boundary. 

The main advantage of the RMT approach over the 
previous implementation is its improved accuracy and 
numerical efficiency. Whereas the previous implemen¬ 
tation used a low-order Crank-Nicolson propagator, the 
RMT approach uses a high-order Arnoldi propagator 
. This replaces a solution of a system of linear equa¬ 
tions by repeated matrix-vector multiplications, which 
may reduce numerical error in the time and spatial prop¬ 
agation of the wavefunction. Since the Arnoldi propa¬ 
gator is dominated by matrix-vector multiplications, the 
RMT codes can be parallelised more efficiently, so that 
calculations exploiting in excess of 2000 cores are feasible. 

In the present study, we aim to investigate inner-shell 
ionization processes involving the carbon atom. We are 
thus interested in residual-ion states with a hole in the 
Is shell. Within the R-matrix codes, residual-ion states 
are retained in order of energy. As a consequence, all 
possible residual-ion states with a filled Is^ shell are in¬ 
cluded prior to inclusion of the residual-ion states with a 
hole in the Is shell. In order to limit the scale of the cal¬ 
culations, we therefore adopt a minimal basis for the de¬ 
scription of carbon. The atom is described using only the 
Is, 2s and 2p Hartree-Fock orbitals of singly ionized car¬ 
bon [ 23 . We then use these orbitals to build all possible 
singly-charged residual-ion states: i.e. all states belong¬ 
ing to the ls^2s^2p, ls^2s2p^, ls^2p^, Is2s^2p^, Is2s2p^ 
and ls2p^ configurations of C’*'. The neutral-atom ba¬ 
sis then contains all combinations of these residual-ion 
states with a set of B-spline based continuum orbitals up 
to a maximum total angular momentum Umax = 5. We 
use a total of 125 B-splines of order 9 to build these con¬ 
tinuum orbitals. This basis also includes so-called corre¬ 
lation functions made of all combinations for 6 electrons 
across the Is, 2s and 2p orbitals with at least one electron 
in Is. The inner-region boundary is set to 27 ag. 

The application of RMT theory using this basis set 
to describe the atom poses a new challenge: a challenge 
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that applies to all studies of inner-shell ionization pro¬ 
cesses. To describe all ionization processes properly, a 
good description of the continuum is needed for emission 
of both inner and outer electrons, including multipho¬ 
ton emission processes. Multiphoton emission of outer 
electrons, in particular, can lead to very high contin¬ 
uum energies. An extensive expansion of the continuum 
is therefore required in the inner region, which allows 
for the description of rapidly oscillating continuum func¬ 
tions. As a consequence, the so-called knot points of the 
B-spline basis set are more closely spaced than in a cal¬ 
culation for outer electrons only: 125 B-splines in the 
present case, compared to 70 for the outer-electron cal¬ 
culation [^. The introduction of the Bloch operator to 
maintain Hermiticity of the Hamiltonian in the inner re¬ 
gion then generates eigenfunctions with large eigenvalues 
(up to 200 keV) which are sharply peaked near the R- 
matrix boundary. Since these functions are peaked near 
the boundary, they need to be retained in the calcula¬ 
tions. However, the effect of the large eigenvalues must 
be compensated for through the connection between the 
inner region and the outer region. Hence, a cancellation 
of terms involving large energies occurs at every stage of 
the calculation. This cancellation can be a prime source 
of numerical error, and extra care therefore needs to be 
taken to ensure numerical stability of the calculations, for 
example through a significant reduction of the time-step 
in the calculation. 

Within the RMT approach, the light field is assumed 
to be linearly polarized and described within the length 
form of the dipole approximation due to the necessity to 
restrict the residual-ion basis [2^ . The field is described 
by an ultra-short light pulse of eight cycles, including a 
three-cycle sin^ ramp-on and ramp-off, with two cycles at 
peak intensity. The photon energies in the present study 
range from 170 eV to 245 eV, so that a single photon 
suffices to eject an outer 21 electron, but absorption of 
two photons is required to emit the Is electron. The 
bandwidth of the pulse is about 40 eV (full width at half 
maximum) at a photon energy of 245 eV. After the pulse 
has ended, we propagate the wavefunction for another 42 
cycles to ensure that all ejected electrons have entered 
the outer region. The time-step used in the calculation 
is 0.012 as. The outer-region finite difference grid has a 
spacing of 0.025 oq, and extends out to a distance of 816 
ao- We use an Arnoldi propagator of order 10. 


III. RESULTS 

In the present study, we aim to investigate the com¬ 
petition between single-photon emission of a 2£ electron 
and two-photon emission of a Is electron from a carbon 
atom in the photon energy range between 170 and 245 
eV. Since the C ground state has symmetry, the ini¬ 
tial state can have = —1,0 and 1. In all results pre¬ 
sented, unless otherwise stated, we have averaged over 
the different initial M^-values. For non-zero initial M^, 


the S symmetry is not available, whereas for zero 
radiative transitions with AL = 0 are not allowed. 

In table [H we present final-state populations in the 
outer region for various subsets of photoemission chan¬ 
nels when carbon is irradiated by a short pulse of 190 
eV and 245 eV photons at various peak laser intensities. 
The table shows that the yield for odd-parity channels 
associated with emission of a 2£ electrons scales approx¬ 
imately linearly with intensity between 10^^ and 10^^ 
W/cm^. This indicates that these channels correspond 
to single-photon emission of an outer 2£ electron. The 
population in the even-parity channels associated with 
2£ emission increases quadratically with intensity, and 
this population can thus be interpreted as two-photon 
above-threshold emission of a 2£ electron. Similarly, the 
population in the even-parity channels associated with Is 
emission increases quadratically with intensity, and this 
population can be interpreted as two-photon emission of 
a Is electron. Up to intensities of 10^^ W/cm“^, the 
assumption that perturbation theory applies is therefore 
not unreasonable, with deviations typically less than 7%. 

At an intensity of 10^® W/cm^, the assumption of per¬ 
turbation theory certainly no longer holds. The single¬ 
photon emission yield for a 2£ electron increases by about 
a factor 6 between 10^^ W/cm^. The two-photon emis¬ 
sion yield for a Is electron increases by about a factor 60 
between 10^^ W/cm^. In both cases, a reduction of about 
40% from the perturbative result is observed. For two- 
photon emission of a 2£ electron, an additional increase 
in the yield of about 10% is seen. Total ionization has be¬ 
come substantial at 10^® W/cm^, with about 27% of all 
atoms ionized within a short time. Hence, perturbation 
should no longer be expected to apply at this intensity. 
Perturbation theory could fail in different respects: an 
increase in importance of higher-order processes, and sat¬ 
uration of ionization, in combination with the very short 
duration of the XUV pulse (smaller than 271/Ip with Ip 
the ionization potential of the 2p electron). Also, the 
strength of the electric held will significantly change the 
actual potential seen by the outer electrons. 

Table U shows, furthermore, that the population in 
odd-parity channels associated with Is emission scales 
linearly with intensity at the lowest intensities for 190 
eV photons and at all intensities for 245 eV photons. 
The reason for this is the large bandwidth of the ultra- 
short light pulse. The threshold for single-photon emis¬ 
sion of the Is electron in the present calculations is about 
300 eV. (This threshold is not necessarily determined 
very accurately, since the lowest Is emission threshold, 
ls2s^2p^ is described using Hartree-Fock orbitals 

for the ls^2s^2p ground state of C’*'.) The ultra-short 
nature of the light pulse has sufficient bandwidth to al¬ 
low single-photon ionization of a Is electron to occur with 
a relative strength, compared to single-photon emission 
of a 2£ electron, of about at 190 eV and about 

3x10“^ at 245 eV. This increase with increasing photon 
energy is expected, as a larger photon energy will lead to 
a larger part of the bandwidth having sufficient energy 
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TABLE I: Final-state populations in the outer region for ground-state C atoms irradiated by an ultra-short laser pulse with 
a central photon energy of 190 and 245 eV at different peak laser intensities. The ground state of C has even parity. The 
notation 1.54(-6) indicates 1.54 x 10“®. The populations are averaged over initial orbital magnetic quantum number Ml- 


Channel subset 

Photon energy 



Peak intensity 




eV 

10^“^ W/cm^ 

10^® W/cm^ 

10^® W/cm^ 

10^^ W/cm^ 

10^® W/cm^ 

C*" 21 emission, odd parity 

190 

3.89(-5) 

3.89(-4) 

3.87(-3) 

3.64(-2) 

2.07(-l) 

C*" 2£ emission, even parity 

190 

2.45(-10) 

2.45^8) 

2.45^6) 

2.45H) 

2.68^2) 

C*" Is emission, even parity 

190 

6.59^10) 

6.59^8) 

6.55^6) 

6.18H) 

3.66^2) 

C"*" Is emission, odd parity 

190 

2.55(-9) 

2.55^8) 

2.57^7) 

7.92^6) 

3.16^3) 

C*" 21 emission, odd parity 

245 

1.29(-5) 

1.28H) 

1.28(-3) 

1.24(-2) 


C'*' 2£ emission, even parity 

245 

5.05(-ll) 

5.05^9) 

5.02^7) 

4.76^5) 


C"*" Is emission, even parity 

245 

3.19^10) 

3.18^8) 

3.16^6) 

2.96H) 


C”'" Is emission, odd parity 

245 

3.90(-7) 

3.90^6) 

3.90^5) 

3.89H) 



to eject a Is electron. At 245 eV, this bandwidth leads 
to dominance of single-photon emission of a Is electron 
over multiphoton emission even at an intensity of 10^® 
W/cm^. However, at 190 eV, the population of the Is 
emission channels with odd parity increases by a factor 
of 400 when the intensity is increased from 10^^ W/cm^ 
to 10^® W/cm^. This indicates that for this photon en¬ 
ergy, three-photon processes start to become important 
near an intensity of 10^^ W/cm^. We can therefore ob¬ 
tain an estimate for three-photon emission of a Is elec¬ 
tron through comparison of the final-state populations 
obtained at intensities of 10^^ W/cm^ and 10^^ W/cm^. 
This comparison gives a final-state population in Is emis¬ 
sion channels associated with three-photon absorption at 
an intensity of 10^^ W/cm^ of around 5.4 x 10“®. This 
procedure to estimate above threshold emission for the 
Is electron can be carried out for photon energies up to 
about 220 eV. 

The final-state populations in table U have been av¬ 
eraged over Ml- Little difference is seen between the 
final-state populations for Ml = 0 and Ml = ±1, except 
for two-photon emission of the 2£ electrons, for which 
the Ml = ±1 yield is about 65% larger than the Ml = 0 
population at 190 eV, a factor 2 larger at 218 eV and 
10% smaller at 245 eV. The most likely reason for the 
generally larger yield for Ml = ±1 is the presence of a 
2 p electron with mt = 0 for the initial Ml = 1 level of 
the ®P® ground state of C, whereas no such electron is 
present for the initial Ml — 0 level. We note that this 
same principle was the reason for the relative increase in 
the emission of 2s electrons compared to 2p electrons in 
multiphoton ionization of C at 390 nm for Ml = 0 0. 

We can compare the ionization yields in the odd-parity 
channels with estimates of the pho toionization cross sec¬ 
tions in the various subshells (29j. The calculation we 
compare with also uses a very simple description for C. 
First we consider the results for a central photon energy 
of 190 eV. The ionization yield for 2£ electrons gives a 
photoionization cross section for the n = 2 subshell of 
0.128 Mb, very similar to the result presented in [^ . 
The determination of a Is photoionization cross section 
is more difficult. The Fourier transform of the electric 
field shows that the photon energy spectrum above the 


Is ionization threshold accounts for about 2.3 x 10“® of 
the full intensity. It contains three main components: 
20 % of the intensity associated with photon energy spec¬ 
trum above the Is ionization threshold is found within 
20 eV of the threshold, 70% of the intensity is associ¬ 
ated with photon energies around 320 eV, and 10% with 
photon energies around 365 eV. If we assume that the 
full Is emission yield arises from photon energies around 
320 eV, a cross section of 0.64 Mb is obtained, which 
compares reasonably well with a reported photoioniza¬ 
tion cross section of 0.86 Mb at a photon energy of 300 
eV [ 2 ^, given the various uncertainties in transforming 
the yield into a cross section. 

Using the same procedure, we obtain a photoioniza¬ 
tion cross section of 0.070 Mb for the n = 2 subshell for 
a central photon energy of 245 eV, in good agreement 
with the result reported in [ 2 ^. The estimate of a Is 
photoionization cross section is again more complicated. 
The Fourier transform of the electric field shows that the 
photon energy spectrum above the Is ionization thresh¬ 
old accounts for 0.88% of the full intensity. It contains 
two main components: 33% of the intensity with photon 
energy above the Is ionization threshold is found within 
8 eV of the threshold, with the remaining 67% of the in¬ 
tensity associated with photon energies peaking around 
310 eV. If we assume that the full Is emission yield arises 
from photon energies around 310 eV, a cross section of 
0.31 Mb is obtained, which is nearly a factor 3 smaller 
than the reported photoionization cross section of 0.86 
Mb at a photon energy of 300 eV [ 2 ^. The reason for 
this increase in discrepancy is unknown. The components 
near threshold are strongest at threshold, and may there¬ 
fore not lead to ionization within the calculations, due to 
the time needed by the photoelectron to leave the inner 
region. However, this would only increase our estimate 
of the cross section to 0.46 Mb, which is still notably be¬ 
low the previous result. Other possible reasons for the 
discrepancy include the short duration of the pulse, and 
possible interference arising from excitation of the Is elec¬ 
tron to np states by the central component of the pulse. 

TableUshows that, through consideration of the outer- 
region channels, we can obtain emission yields for single¬ 
photon emission of a 2s or 2p electron, two-photon emis- 





5 



FIG. 1: (Colour online) Yields for one- and two-photon emis¬ 
sion of a electron and two- and three-photon emission of 
a Is electron from neutral C irradiated by a short laser pulse 
with a peak intensity of W/cm^ as a function of photon 
energy. The three-photon emission yield is estimated through 
subtraction of the single-photon yield induced by the edges of 
the pulse bandwidth. Data associated with this report can be 
accessed via [^ . 

sion of a Is electron, two-photon emission of a 2s or 2p 
electron, and an approximation to three-photon emission 
yield of a Is electron at lower photon energies at an in¬ 
tensity of 10^’’ W/cm^. Figure [T] shows these emission 
yields as a function of photon energy. The figure demon¬ 
strates that multiphoton emission is not negligible at this 
intensity, as it contributes between 2 and 3% to the total 
emission. The figure demonstrates that two-photon emis¬ 
sion of the inner Is electron is more important than two- 
photon emission of an outer 2s or 2p electron, and the im¬ 
portance of two-photon emission of a Is electron increases 
with photon energy, compared to both single-photon and 
two-photon emission of a 2^ electron. Three-photon emis¬ 
sion of the Is electron contributes about 0.02% to the 
total photoemission yield for a photon energy of 170 eV 
and an intensity 10^^ W/cm^. This contribution drops 
to about 0.01% at a photon energy of 210 eV. 

At the highest intensity considered in the present cal¬ 
culation, 10^® W/cm^, the increase in ionization yield no 
longer follows a perturbative scaling law. However, the 
relative importance of higher-order processes continues to 
increase. At this intensity, two-photon emission of a Is 
electron accounts for 16% of all photoemission processes, 
with two-photon emission of a 2£ electron accounting for 
another 2%. Therefore at intensities approaching 10^® 
W/cm^, a full analysis of the physics must account for 
the potential effects of absorption of more than just a 
single photon. 

Previously, we investigated the competition between 
multiphoton absorption by a Is and 2s electron from the 
initial ls2s state in He [l^. In that study, it was 
found that at the two-photon level, the Is electron was 
about five times as likely as the 2s electron to absorb 
two photons. In the present study, at 175 eV, the Is 
electron is about 2.5 times as likely to absorb two photons 


TABLE II: Populations of final states of after irradiation 
of a ground-state C atom by an eight-cycle laser pulse with 
an intensity of 10^^ W/cm^ at photon energies of 170 eV, 197 
eV and 224 eV. The notation 1.16(-2) denotes 1.16 x 10~^. 


c+ 

Final-state population 


at a photon energy of 

state 

170 eV 

197 eV 

224 eV 

T?2?2p^P' 

1.16(-2) 

5.35(-3) 

2.64(-3) 

ls''2s2p^ ""P® 

3.07(-2) 

1.70(-2) 

9.98(-3) 

Is^2s2p2 

4.92(-4) 

2.12(-4) 

l.ll(-4) 

Is22s2p2 2s« 

7.47^5) 

3.43^5) 

2.03(-5) 

Is^2s2p2 2p® 

1.24(-2) 

7.24(-3) 

4.43(-3) 

ls^2p'' '‘8° 

9.92(-4) 

5.07^4) 

2.88(-4) 


1.55^3) 

7.81^4) 

4.29^4) 

Is2 2p3 2po 

8.29H) 

4.40^4) 

2.51^4) 

ls2s^2p^ ""P® 

4.90(-4) 

3.56(-4) 

2.29(-4) 

Is2s^2p2 2p® 

2.91H) 

1.74(-4) 

1.08^4) 

Is2s^2p2 

9.27^6) 

8.80^7) 

2.52(-6) 

Is2s^2p2 

3.18^6) 

1.61^7) 

3.58^7) 

Is2s2p'' -"8° 

2.71^6) 

1.21(-6) 

5.01(-6) 

Is2s2p® "*0° 

9.94(-6) 

8.06^6) 

8.43^6) 

Is2s2p® '‘P° 

4.47^6) 

3.60^6) 

4.60^6) 

Is2s2p® ^D° 

1.12(-5) 

4.66^6) 

2.85^6) 

Is2s2p® ^P° 

5.46(-6) 

1.85^6) 

1.30^6) 

Is2s2p® ^D° 

3.51^6) 

1.03^6) 

2.63^6) 

Is2s2p® ■‘S" 

4.29(-6) 

3.03^6) 

1.80(-6) 

Is2s2p® ^S° 

1.62(-6) 

1.23^6) 

2.48^6) 

Is2s2p® ^P° 

3.11^6) 

5.17^7) 

1.49^6) 

ls2p^ ^P® 

4.59^7) 

2.85^7) 

3.37^7) 

ls2p'‘ 

7.30^7) 

1.79^7) 

3.48^7) 

ls2p'‘ 

3.48^7) 

7.21^8) 

2.09^7) 

ls2p'‘ 2p® 

2.11(-7) 

3.17^8) 

3.49^8) 


than for a ‘It electron to absorb two photons, whereas it 
is a factor 6 at a photon energy of 245 eV. The ratio 
obtained here is therefore comparable to the one found 
for He, even though in the present case, there are twice as 
many 2£ electrons as Is electrons. On the other hand, as 
shown later, the emission of 21 electrons is dominated by 
emission of a 2s electron, with 2p electrons contributing 
notably less. 

The RMT calculations also provide the populations in 
the various residual 0+ states. Table [III provides the final 
populations for all 0+ states included in the calculations 
for photon energies of 170 eV, 197 eV and 224 eV at an 
intensity of 10^^ W/cm^. For the photon energies given 
in the table, effects due to the bandwidth of the pulse are 
expected to play only a minor role, and should not affect 
the final populations significantly. 

The table demonstrates that the main residual ion 
states are the ls^2s2p^ and states indicating a 
dominance of direct emission of a 2s electron. Emission 
of a 2p electron accounts for 20% of the emission pro¬ 
cesses at 170 eV and only 14% at 224 eV. C+ is left in a 
ls^2p® state in about 5.2-5.7% of all emission processes. 
This latter outcome is only a factor 4 more likely than the 
emission of an inner Is electron, which occurs in 1.4-2.0% 
of all emission processes. 
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FIG. 2: (Colour online) Photoionization yields for residual- 
ion configurations of as a function of photon energy at an 
intensity of 10^^ W/cm^. 

In addition, table HIl shows that the emission of the in¬ 
ner electrons can involve the outer electrons. The outer 
electrons are left in a 2s^2p^ configuration in 91 - 93% of 
all Is emission processes. However, the final distribution 
over the ls2s^2p^ shows deviations from the statistical 
distribution of a 2:1 ratio between the and the 
state, with a 5:3 ratio at 170 eV and a ratio of 2.1:1 at 
a photon energy of 224 eV. In 7-9% of the Is ionization 
processes, a change in the outer electron population oc¬ 
curs, mainly an excitation of a 2s electron to 2p. This 
suggests that in just under 10% of the processes, two dif¬ 
ferent electrons absorb a photon. Hence this provides a 
signature of a multi-electron response to the light field, 
which may involve a sequential process, whereby the first 
photon excites a Is electron to 2p, followed by photoe¬ 
mission of one of the 2s electrons. 

The table demonstrates that all emission processes in¬ 
volving an outer 2s or 2p electron decrease in magni¬ 
tude with increasing photon energy. However, this pat¬ 
tern changes when a Is electron is emitted. Although 
the population in the dominant residual-ion states after 
emission of a Is electron (ls2s^2p^ ^P® and '^P^) shows 
a decrease in final population with increasing photon en¬ 
ergy, for the other ls2s^2p^ states, for five out of nine 
Is2s2p^ residual-ion states and all ls2p‘* states, the final 
population is larger at a photon energy of 224 eV than 
at a photon energy of 197 eV. This could be due to an 
intermediate resonance, but it could also originate from 
the closer proximity of the threshold for single-photon 
emission of a Is electron at this higher photon energy. 

To investigate the origin of the increase, we have car¬ 
ried out further calculations of inner-shell photoemission 
of C atoms covering the photon-energy range between 170 
and 245 eV. Outcomes of these calculations are shown in 
figured To maintain clarity, the figure shows the proba¬ 
bility that the 0“'" ion is left in a particular configuration, 
rather than in individual states within the configuration. 


Most residual O’*" ions are left in the ls^2s2p^ configura¬ 
tion, indicating that emission of a 2s electron is the most 
likely photoionization process. The probability of emis¬ 
sion of a 2p electron is about a factor 4 times smaller 
at 170 eV, and this probability decreases with increas¬ 
ing photon energy. The probability that the 0“'" ion is 
left in ls^2p^ is over an order of magnitude smaller than 
the probability for emission of a 2s electron. Emission 
of a Is electron preferentially leaves the C’*' ion in the 
ls2s^2p^ configuration, accounting for about 2% of the 
photoionization processes at 170 eV. This probability in¬ 
creases beyond 230 eV, and reaches about 5% for 245 
eV. This increase at the highest photon energies can be 
ascribed to the increasing proximity of the threshold for 
single-photon emission of the Is electron. 

Figure[2]shows that the probability of leaving the resid¬ 
ual C’*' ion in the Is2s2p^ and ls2p'* conhgurations are 
resonantly enhanced at a photon energy of 224 eV. The 
populations in the ls2s^2p^ configuration show an en¬ 
hancement near the threshold for single-photon emission 
of a Is electron, but this enhancement only becomes no¬ 
ticeable above a photon energy of 230 eV. The most likely 
origin of the enhancement is the intermediate ls2s^2p^ 
configuration reached by photoexcitation of a Is electron 
into the 2p shell. The ultra-short duration of the pulse 
means that the effect of the intermediate configuration 
is apparent over a broad photon energy range with very 
little structure. Since the effects of the resonances are 
spread out over a broad photon energy range, the overall 
increase in the ionization probability will be reduced, and 
as a consequence resonant enhancement is only visible for 
the weakest photoionization channels. 


IV. CONCLUSIONS 

In conclusion, we have demonstrated the capability of 
RMT theory to investigate ultra-fast inner-shell emission 
processes in general multi-electron atoms. Two-photon 
emission of a Is electron from C atoms was investigated 
in the photon energy range between 175 and 245 eV. At 
an intensity of 10^^ W/cm^, two-photon emission of the 
Is electron accounts for about 2-3% of all photoionization 
processes. At the two-photon level, emission of the Is 
electron is about a factor two more likely than the emis¬ 
sion of either a 2p or a 2s electron combined. Through 
examination of the final configuration of the residual C“'" 
ion, we have furthermore observed evidence for resonant 
enhancement of the Is emission. Finally, by comparing 
ionization yields in different symmetries across multiple 
intensities, we have determined an above-threshold ion¬ 
ization yield for inner-shell photoemission of the Is elec¬ 
tron in C. 

The determination of inner-shell processes in general 
multi-electron atoms poses a number of challenges. The 
present calculations have been carried out using Hartree- 
Fock orbitals for ground-state C. The description of 
the residual ion states can be improved by expanding 
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the initial orbital set to include pseudo orbitals. How¬ 
ever, the inclusion of pseudo-orbitals will lead to signif¬ 
icantly larger calculations. The main reason for this is 
that residual-ion states are retained in the present R- 
matrix calculations in order of energy. The lowest en¬ 
ergy states will be dominated by states with two Is elec¬ 
trons. These states can involve the Hartree-Fock 2s and 
2 p orbitals, but may also involve the additional pseudo¬ 
orbitals. Hence, a significant number of states dominated 
by pseudo-orbitals may have to be included before one 
reaches the ionic states corresponding to emission of a Is 
electron. Furthermore, the calculations require a good 
description of the continuum up to high energy. This 
leads to much higher energies associated with the Bloch 
operator. The higher energies reached in the present cal¬ 
culations do not yet affect the stability of the RMT ap¬ 
proach. 

A further open question concerns the influence of non¬ 
dipole terms. The present study assumes that the dipole 


approximation holds. This approximation is not unrea¬ 
sonable for the interaction between the laser field and 
the Is orbital. However, the approximation may not be 
as suitable for describing the interactions between the 2£ 
electrons and the laser field. The investigation of non¬ 
dipole effects would require substantial computational 
development. The RMT codes themselves will need to 
be modified to ensure that all relevant interactions aris¬ 
ing from the laser field are properly taken into account. 
More importantly, the inner-region R-matrix codes need 
modification so that higher-order transitions are calcu¬ 
lated within the inner region, and generated in a form 
suitable for use within the RMT codes. 
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